
Nonlinear Optics 5A5513 (2003) 
Lecture notes 



Lecture 7 



So far, this course has mainly dealt with the dependence of the angular frequency of the light 
and molecular interaction strength in the description of nonlinear optics. In this lecture, we will 
now end this development of the description of interaction between light and matter, in favour 
of more engineering practical techniques for describing the theory of an experimental setup in a 
certain geometry, and for reducing the number of necessary tensor elements needed for describing 
a mediumof a certain crystallographic point-symmetry group. 

Motivation for analysis of susceptibilities in rotated coordinate systems 

For a given experimental setup, it is often convenient to introduce some kind of reference coordinate 
frame, in which one for example express the wave propagation as a linear motion along some 
Cartesian coordinate axis. This laboratory reference frame might be chosen, for example, with 
the z-axis coinciding with the direction of propagation of the optical wave at the laser output, in 
the phase-matched direction of an optical parametric oscillator (OPO), after some beam aligning 
mirror, etc. 

In some cases, it might be so that this laboratory frame coincide with the natural coordinate 
frame^ of the nonlinear crystal, in which case the coordinate indices of the linear as well as nonlinear 
susceptibility tensors take the same values as the coordinates of the laboratory frame. However, we 
cannot generally assume the coordinate frame of the crystal to coincide with a conveniently chosen 
laboratory reference frame, and this implies that we generally should be prepared to spatially 
transform the susceptibility tensors to arbitrarily rotated coordinate frames. 

Having formulated these spatial transformation rules, we will also directly benefit in another 
aspect of the description of nonlinear optical interactions, namely the reduction of the suscepti- 
bility tensors to the minimal set of nonzero elements. This is typically performed by using the 
knowledge of the so called crystallographic point ymmctry group of the medium, which essentially 
is a description of the spatial operations (rotations, inversions etc.) that define the symmetry 
operations of the medium. 

As a particular example of the applicability of the spatial transformation rules (which we 
soon will formulate) is illustrated in Figs. 1 and 2. In Fig. 1, the procedure for analysis of sum 
or frequency difference generation is outlined. Starting from the description of the linear and 
nonlinear susceptibility tensors of the medium, as we previously have derived the relations from 
a first principle approach in Lectures 1-6, we obtain the expressions for the electric polarization 
densities of the medium as functions of the applied electric fields of the optical wave inside the 
nonlinear crystal. These polarization densities arc then inserted into the wave equation, which 
basically is derived from Maxwell's equations of motion for the electromagnetic field. In the wave 
equation, the polarization densities act as source terms in an otherwise homogeneous equation for 
the motion of the electromagnetic field in vacuum. 

As the wave equation is solved for the electric field, here taken in complex notation, we have 
solved for the general output from the crystal, and we can then design the experiment in such a 
way that an optimal efficiency is obtained. 



^ The natural coordinate frame of the crystal is often chosen such that some particular symmetry 
axis is chosen as one of the Cartesian axes. 
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Laboratory frame (a;', j/', z') 
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(Nonlinear wave equation inside crystal) 



E„. 



/c = 1,2,3,4 



(Solution for optical fields) 
Figure 1. The setup in which the orientation of the laboratory and crystal frames coincide. 

In Fig. 1, this outline is illustrated for the case where the natural coordinate frame of the crystal 
happens to coincide with the coordinate system of the laboratory frame. In this case, all ele- 
ments of the susceptibilities taken in the coordinate frame of the crystal (which naturally is the 
coordinate frame in which we can obtain tabulated sets of tensor elements) will coincide with the 
elements as taken in the laboratory frame, and the design and interpretation of the experiment is 
straightforward. 

However, this setup clearly constitutes a rare case, since we have infinitely many other possibil- 
ities of orienting the crystal relative the laboratory coordinate frame. Sometimes the experiment is 
designed with the crystal and laboratory frames coinciding, in order to simplify the interpretation 
of an experiment, and sometimes it is instead necessary to rotate the crytal, in order to achieve 
phase-matching of nonlinear process, as is the case in for example most schemes for second-order 
optical parametric amplification. 

If now the crystal frame is rotated with respect to the laboratory frame, as shown in Fig. 2, we 
should make up our mind in which system we would like the wave propagation to be analyzed. In 
some cases, it might be so that the output of the experimental setup is most easily interpreted in 
the coordinate frame of the crystal, but in most cases, we have a fixed laboratory frame (fixed by 
the orientation of the laser, positions of mirrors, etc.) in which we would like to express the wave 
propagation and interaction between light and matter. 
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Laboratory frame {x',y',z') Crystal frame {x,y,z) [rotated relative {x',y',z')] 



^ Ei(r, t) = Y:1=i Re[E„, exp{-iujkt)] ^ 



huJ2 

WW- 



Al) , (1) (2) , (2) 
X.x'x' r Xxx , Xx'y'z' r Xxyz, 

(Properties of the medium) 




-Moa;2(PLV+PL'2)=? fc = 1,2,3,4 



(Nonlinear wave equation inside crystal) 



(Solution for optical fields) 

Figure 2. The setup in which the crystal frame is rotated relative the laboratory frame. 
In Fig. 2, wc would, in order to express the nonlinear process in the laboratory frame, like to 

(2) (2) 

obtain the naturally appearing susceptibilities Xxyz, Xxxx, etc., in the laboratory frame instead, as 
xi'l/'z', xi'i'x', etc. 

Just to summarize, why are then the transformation rules and spatial symmetries of the meduim 
so important? 

• Hard to make physical conclusions about generated optical fields unless orientation of the 
laboratory and crystal frames coincide. 

• Spatial symmetries often significantly simplifies the wave propagation problem (by choosing a 
suitable polarization state and direction of propagation of the light, etc.). 

• Useful for reducing the number of necessary elements of the susceptibility tensors (using Neu- 
mann's principle). 

Optical properties in rotated coordinate frames 

Consider two coordinate systems described by Cartesian coordinates Xa and x'^, respectively. The 
coordinate systems are rotated with respect to each other, and the relation between the coordinates 
are described by the [3 x 3] transformation matrix Rat as 

x' = Rx ^ x'„ = Rafjxp, [B. & C. (5.40)] 

where x = (x, y, z)^ and x' = (a;', y' , z')^ are column vectors. The inverse transformation between 
the coordinate systems is similarly given as 

x = R-ix' <^ Xf) = R„0x'^. [B.&C. (5.41)] 
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Figure 3. Illustration of proper rotation of the crystal frame [x, y, z) relative to the laboratory 
reference frame {x', y', z'), by means of x'^ = RapXp with det (R) = 1. 




(b) 

z 




z' 



Figure 4. The coordinate transformations (a) x = {x,y,z) hh. x' = (— a;,— y, z), constituting a 
proper rotation around the z-axis, and (b) the space inversion x i— » x' = — x, an improper rotation 
corresponding to, for example, a rotation around the z-axis followed by an inversion in the xy-plane. 

We should notice that there are two types of rotations that are encountered as transformations: 

• Proper rotations, for which det(R) = 1. (Righthanded systems keep being righthanded, and 
lefthanded systems keep being lefthanded.) 

• Improper rotations, for which det(R) = — 1. (Righthanded systems are transformed into left- 
handed systems, and vice versa.) 

The electric field E(r, €) and electric polarization density P(r, i) are both polar quantities that 

transform in the same way as regular Cartesian coordinates, and hence we have descriptions of 
these quantities in coordinate systems {x,y,z) and {x',y',z') related to each other as 

E'^ir, t) = R^uEuir, t) ^ E^{r, t) = R^^uE'^ir, t), 

and 

P;(r,t) = R^uPu{r,t) ^ Pu{r,t) = i?^„P;(r,t), 

respectively. Using these transformation rules, we will now derive the form of the susceptibilities 
in rotated coordinate frames. 

First order polarization density in rotated coordinate frames 

From the transformation rule for the electric polarization density above, using the standard form as 
we previously have expressed the electric field dependence, we have for the first order polarization 
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density in the primed coordinate system 

P«'(r,t) = i?^„P«(r,t) 

/OO 
-OO 

/OO 
xia (-w; a;)-Raa-E'a (w) exp(-ia;t) cL; 
-OO 

/>oo 
-OO 

where 

xW'(-a;; a;) = i?^„i?„„xl'i (-c^; [B. & C. (5.45)] 

is the Hnear electric susceptibility taken in the primed coordinate system. 
Second order polarization density in rotated coordinate frames 

Similarly, we have the second order polarization density in the primed coordinate system as 
Pf (r,i) = i?^„p(2)(r,t) 

/■OO /-CO 
/ Xulbi~'^<^'i^i^'^2)Ea{oJi)Eb{uj2)exp[-i{oJi + u>i)t] dijj2 dui 
-OC J — OO 

/OO />oo 
/ xia6(-'^<^;'^i''^2)-Rao-E^('^i)^/36-E^(w2)exp[-i(a;i +a;i)t]rfw2rfa;i 
-OO J — OO 

/OO pOO 
/ x[,a/3'(-^cr; wi, W2)-E^(wi)-E^(a;2) exp[-i(a;i + u)i)t] du}2 dwi 
-OO J —OO 

where 

x|ta/3'(-Wa;Wl,a;2) = i?^„-Raa-R/35xiafe(-'^<^;'^l'^2) [B. & C. (5.46)] 

is the second order electric susceptibility taken in the primed coordinate system. 

Higher order polarization densities in rotated coordinate frames 

In a manner completely analogous to the second order susceptibility, the transformation rule be- 
tween the primed and unprimed coordinate systems can be obtained for the nth order elements of 
the electric susceptibility tensor as 

{-u;y,ui,...,iVn). [B.&C. (5.47)] 



Crystallographic point symmetry groups 

Typically, a particular point symmetry group of the medium can be described by the generating 
matrices that describe the minimal set of transformation matrices (describing a set of symmetry 
operations) that will be necessary for the reduction of the constitutive tensors. Two systems are 
widely used for the description of point symmetry groups:^ 

• The International system, e. g. 43m, m3m, 422, etc. 

• The Schonflies system, e. g. T^, Oh, D4, etc. 

The crystallographic point symmetry groups may contain any of the following symmetry oper- 
ations: 

1. Rotations through integral multiples o/27r/n about some axis. The axis is called the n-fold 
rotation axis. It is in solid state physics shown [1-3] that a Bravais lattice can contain only 2-, 

3-, 4-, or 6-fold axes, and since the crystallographic point symmetry groups are contained in the 
Bravais lattice point groups, they too can only have these axes. 

^ C. f. Table 2 of the handed out Hartmann's An Introduction to Crystal Physics. 
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2. Rotation-reflections. Even when a rotation through 27r/n is not a symmetry element, 
sometimes such a rotation followed by a reflection in a plane perpendicular to the axis may be a 
symmetry operation. The axis is then called an n-fold rotation-reflection axis. For example, the 
groups 5*6 and S'4 have 6- and 4-fold rotation-reflection axes. 

3. Rotation-inversions. Similarly, sometimes a rotation through 27r/n followed by an inversion 
in a point lying on the rotation axis is a symmetry element, even though such a rotation by itself 
is not. The axis is then called an n-fold rotation-inversion axis. However, the axis in Sq is only a 
3-fold rotation-inversion axis. 

4. Reflections. A reflection takes every point into its mirror image in a plane, known as a 
mirror plane. 

5. Inversions. An inversion has a single fixed point. If that point is taken as the origin, then 
every other point r is taken into — r. 

Schonflies notation for the non-cubic crystallographic point groups 

The twenty-seven non-cubic crystallographic point symmetry groups may contain any of the fol- 
lowing symmetry operations, here given in Schonflies notation^: 

Cn These groups contain only an n-fold rotation axis. 

C„v In addition to the n-fold rotation axis, these groups have a mirror plane that contains 
the axis of rotation, plus as many additional mirror planes as the existence of the 
n-fold axis requires. 

Cnh These groups contain in addition to the n-fold rotation axis a single mirror plane 
that is perpendicular to the axis. 

Sn These groups contain only an n-fold rotation-reflection axis. 

Dn In addition to the n-fold rotation axis, these groups contain a 2-fold axis perpendic- 
ular to the n-fold rotation axis, plus as many additional 2-fold axes as are required 
by the existence of the n-fold axis. 

Dnh These (the most symmetric groups) contain all the elements of Z)„ plus a mirror 
plane perpendicular to the n-fold axis. 

Dnd These contain the elements of D„ plus mirror planes containing the n-fold axis, which 
bisect the angles between the 2-fold axes. 

Neumann's principle 

Neumann's principle simply states that any type of symmetry which is exhibited by the point 
symmetry group of the medium is also possessed by every physical property of the medium. 

In other words, we can reformulate this for the optical properties as: the susceptibility tensors 
of the medium must be left invariant under any transformation that also is a point symmetry 
operation of the medium, or 

(n)' /• (n) / N 

X([*ai--a„(,— <^ct; Wi, . . . ,(J„) — Xnai — a„\—(^a',^l, ■ ■ ■ ,l^n), 

where the tensor elements in the primed coordinate system are transformed according to 

Xfiai-.-an ( '^CT) ^1) • • • ) ^n) ~ ^^u-^aiai ' ' ' ^andnXuai-.-an ( ^a'l ^1 ? * • * ? '^n)? 

where the [3 x 3] matrix R describes a point symmetry operation of the system. 



^ In Schonflies notation, C stands for "cyclic", D for "dihedral", and S for "spiegel". The 
subscripts h, v, and d stand for "horizontal", "vertical", and "diagonal", respectively, and refer 
to the placement of the placement of the mirror planes with respect to the n-fold axis, always 
considered to be vertical. (The "diagonal" planes in Dnd are vertical and bisect the angles between 
the 2-fold axes) 
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Inversion properties 

If the coordinate inversion Ra/j = —Sap , is a symmetry operation of the medium (i. e. if the medium 
possess so-called inversion symmetry), then it turns out that 

(n) n 

for all even numbers n. (Question: Is this symmetry operation a proper or an improper rotation?) 
Euler angles 

As a convenient way of expressing the matrix of proper rotations, one may use the EuJer angles of 
classical mechanics,^ 

where 

(cost/' sint/iOX /lO ^ \ ( cos^ sin(p 0\ 

— sin^ cosV' I , B(t?) = I cosi? sini? I , C(<p) = I —smip cosip I . 

1/ \0 -sini? cost?/ \ 1/ 

Example of the direct inspection technique applied to tetragonal media 

Neumann's principle is a highly useful technique, with applications in a wide range of disciplines in 
physics. In order to illustrate this, we will now apply Neumann's principle to a particular problem, 
namely the reduction of the number of elements of the second order electric susceptibility tensor, 
in a tetragonal medium belonging to point symmetry group 422. 




422 



Figure 5. An object^ possessing the symmetries of point symmetry group 422. 

By inspecting Tables 2 and 3 of Hartmann's An introduction to Crystal Physics^ one find that the 
point symmetry group 422 of tetragonal media is described by the generating matrices 



and 




twofold rotation 
about xi axis 



fourfold rotation 
about 0:3 axis 



C. f. Herbert Goldstein, Classical Mechanics (Addison- Wesley, London, 1980). 

^ The figure illustrating the point symmetry group 422 is taken from N. W. Ashcroft and 
N. D. Mermin, Solid state physics (Saunders College Publishing, Orlando, 1976), page 122. 

^ Ervin Hartmann, An Introduction to Crystal Physics (University of Cardiff Press, International 
Union of Crystallography, 1984), ISBN 0-906449-72-3. Notice that there is a printing error in 
Table 3, where the twofold rotation about the Xs-axis should be described by a matrix denoted 
"M2", and not "Mi" as written in the table. 
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Does the 422 point symmetry group possess inversion symmetry? 

In Fig. 6, the steps involved for transformation of the object into an inverted coordinate frame are 
shown. 



X3 X3 




Figure 6. Transformation into an inverted coordinate system (a;", y" , z") = {—x, —y, —z). 

The result of the sequence in Fig. 6 is an object which cannot be reoriented in such a way that one 

obtains the same shape as we started with for the non-inverted coordinate system, and hence the 
object of point symmetry group 422 does not possess inversion symmetry. 



Step one - Point symmetry under twofold rotation around the xi-axis 

Considering the point symmetry imposed by the R = M4 matrix, we find that (for simplicity 
omitting the frequency arguments of the susceptibility tensor) the second order susceptibility in 
the rotated coordinate frame is described by the diagonal elements 



3 3 3 

= X! X! Rll^RlaRiax'f^ap 

fx=l a = l 13=1 

3 3 3 

= ^ii^^ic.Si0X% = Xm, (identity) 

^=la=l/3=l 



and 

X222 = -R2M-R2Q-R2/3X^a/3 

3 3 3 

= Z Z Z ^2/*^2a-R2/3x|fc^/3 

= E E E(-'52M)(-^2a)(-<52/3)xS^ = -X^^2 
(n=la=l/3=l 

= {Neumann's principle} = X222 = 

which, by noticing that the similar form R^a = ^S^a holds for the 333-component (i. e. the 
2;2;2:-component), also gives xsss = — Xsss = 0- Further we have for the 231-component 



X23I = R2f,R3aRipX^^^l(3 
3 3 3 

= E E E R■2^lR3aRl|3X^^l|} 

^=1 a=l 13=1 

3 3 3 

= E E E(~^2,^)(~^3«)^i/3^i/5 = ^^^ssi, (identity) 

^=la=l /3=1 
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etc., and by continuing in this manner for all 27 elements of X^^^^, one finds that the symmetry 
operation R = M4 leaves us with the tensor elements listed in T?able 1. 



Zero elements Identities (no further info) 

(2) (2) (2) (2) 
Xll2' X1131 X12I1 X13I' 

X2ii> X222> X223) X232 ) (aU Other 13 elements) 

(2) (2) (2) (2) 
X233' X3II' X322' X323> 
(2) (2) 
X332 5 X333 



Table 1. Reduced set of tensor elements after the symmetry operation R = M4. 
Step two ~ Point symmetry under fourfold rotation around the x^-axis 

Proceeding with the next point symmetry operation, described by R = M7, one finds for the 
remaining 13 elements that, for example, for the 123-element 

Xm = RltiR2aRif3X^i2tD 
3 3 3 

= EEE^1^^2ai?3/3xS/3 

^=1 a=l 13=1 

3 3 3 

= E E E(-'^2M)^la^3/3x2/3 = -X213 
^=1 c«=l 13=1 

= {Neumann's principle} = ~Xa2Zi 

and for the 132-element 

Xl32 = RlnR3aR2l3X^i^lp 
3 3 3 

= EEE^1^^3ai?20x!fi 

fj,= l a=l (3=1 

3 3 3 

= E E E(-'^2^)^3a^l/3xS/3 = -X23I 

^=1 a=l 13=1 

= {Neumann's principle} = xisl' 

while the lll-clcmcnt (which previously, by using the R = M4 point symmetry, just gave an 
identity with no further information) now gives 

Xm = -RiM^ia-Ri/3x2/3 

3 3 3 

= E E E -'^l/'^l«-^l/3>^!fa/3 
iu=l c«=l /3=1 

= E E E(-'52^)(-'52a)(-.52^)xS^ = -X'il 
lj.=l a=l /3=1 

= {from previous result for X222} = 
= {Neumann's principle} = Xm • 

By (again) proceeding for all 27 elements of x^a/j) one finds the set of tensor elements as listed in 
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Table 2. (See also the tabulated set in Butcher and Cotter's book, Table A3. 2, page 299.) 



Zero elements Nonzero elements 



(2) (2) (2) (2) (2) (2) _ _ (2) 

Xlll; Xi12j Xi13' Xi21' Xl22' Xl23 ~ X213' 

(2) (2) (2) (2) (2) (2) _ _ (2) 

Xiai) Xi33j X2111 X212' X221) X132 — X231' 

(2) (2) (2) (2) (2) (2) _ _ (2) 

X222> X223' X232' X233' X31IJ X32I — X3I21 

(2) (2) (2) (2) (2) (2) ,„ o • J J 

X313. X322. X323. X331. X332' X333 (6 nonzero, 3 mdependent) 



Table 2. Reduced set of tensor elements after symmetry operations R = M4 and R = M7. 
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